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ABSTRACT 


The influence of internal heating by mesoscale eddies on the 
Oceanic temperature structure is tested by including a vertical heat 
transport term, w't', in the ocean thermal budget equation. The 
value for this term was determined from a theoretical study by Gill, 
Green and Simmons (1974). In the past studies of the thermocline 
problem, the heat transport term has been neglected. Analytical and 
numerical test cases were run solving the heat equation with and 
without the mesoscale eddy heat transport term for baroclinic and 
barotropic conditions in the ocean. Whenever heating by mesoscale 
eddies was included in the heat equation, excessive cooling occurred 
in the upper 1000 m of the ocean. The results of the test cases in- 
dicate that the eddies do have a distinct influence on the ocean's 
thermal structure and need to be included in solutions to the heat 


equation for the ocean. 


TABLE OF CONTENTS 


Tee MIUERODUCTION = - = - — === - = 55 ----=------ 9 
II. EQUATION AND MODEL DEVELOPMENT - - ------------ i) 
A. HEAT EQUATION- - ---------+----+------ 13 

B. M EQUATION ---------------------- 14 

C. BOUNDARY CONDITIONS- - - rf fefrrrererrerrrtcco 16 

Pee oy ee een Se = = = 18 

III. ANALYTIC AND NUMERICAL RESULTS - -----+--+-e+e2+f+-- 2 23 
A. VERTICAL MOTION AT THE SURFACE EQUAL TO ZERO - ---- 23 

l. Analytical Solution- - -----+-+e-fe8-+ef%e- 23 

2. Numerical Solution ----+-+-+-+f+ +++ 72°--- 24 

B. INTERNAL FORCING EQUAL TO ZERO ------------ 27 

1. Analytical Solution- ------+e*4-e8-+f+f+e 2% ZY 

2. Numerical Solution - --------+ +++ 722 30 

C. GENERAL CASE -----+-+-e+-e+-f+e+++7277-°72---°- 33 

D. BAROCLINIC CASES - -----+-+--+2+-+-272+e+2+2+ee-6 a 

1. Internal Heating Equal to Zero -----+-+-+--6- 35 

a. Positive Temperature Gradient- - ------- 38 

b. Negative Temperature Gradient- - ----+-+-- 38 

2. Influence of Internal Heating- - --------- 43 

IV. CONCLUSIONS- - -- -f-r fff frrrrrrrrcrrfrr rece 50 
LIST OF REFERENCES - ---+cfctfffrerr fr ff fer cer fr fr eco 52 
INITIAL DISTRIBUTION LIST- - -]|- -=frfrefrrererrrrrrreree 53 





LOE 


ll. 


1A 


LIST OF FIGURES 


Two dimensional model along 36°N with Az = 40m and 
Ax = 80km ------------------------- ey 


Profile of the rate of conversion of available 
potential energy to kinetic energy. From Gill 
Green and Simmons (1974)- ----+-e+e-eeef8--8--++%+2-- 20 


Mean velocity profiles from Gill, Green and Simmons 
(1974) used to determine temperature profiles - ----- - au 


Profile of vertical derivative of w'T' showing a 
large upward flux in the first 500m. Units are 
°c sect x 1077 -- -------------------- 22 


Analytic solution to the case in which wind forcing 
is equal to zero- ----------f--f- ee ee er eee oe 25 


Initial temperature (°C) profile for numerical model- - - - 26 


Intermediate numerical solution to case in which 
wind forcing equal to zero. a) 28 years pb) 55 
years c) 82 years d) 165 years- -------+------ 28 


Analytic solution to the case in which heating by meso- 
scale eddies is equal to zero with w= 10 cm day"! - - - - 31 


Note the O° lines (no temperature gradient) which 

represent the barotropic solution to the balance 

between w(0T/dz) and K(0°T/dz*). From Stommel 

and Webster (1962)--------+ +--+ e-+e-° 32 


Numerical solution to the case in which heating 
by mesoscale eddies is equal to zero with 
ae -10 cm day71 -~---- ee eee eer er ere ere ee ee 34 


Numerical solution to the general case with 
wy, = 10 cm @day7l --------------------- 36 


Numerical solution to the general case with 
w, = -10 cm day" -------+---+---+----+--+--+---- ay 


We 





Lee 


14. 


LSA 


16. 


lager 


Si 


Ie )e 


ZO. 


On 


Initial temperature field with a 5°C positive 


gradient in the eaSt-west direction- - -------- 


Numerical solution for (dT/dx), > O and lies 


10 cm day +, (0/dz) (w'T') = O0------+------- 


Numerically calculated internal vertical motion 
mor (dT/ox) > O and aes 10 cm day~1. wis in 


cm sec7l x TOs? SoS 6 Fe SS SS SS SS SS 


Initial temperature field with a 5°C negative 


gradient in the east-west direction- - -------- 


Numerically calculated internal vertical motion 
in cm sec7+ x 107° for (OT/0x), < O and Wo = 


10 cm day~1 Se ate ae = Se Ste ee ene eS See Se 


Internal vertical motion for (0T/0x),< O and 


oe -10 cm day~1, w is in cm sec™+ x 107°°- ----- 


Numerical solution for (dT/dx)_< O and wo <A02 


Temperature is in °C. (d/dz) (w'T')= O-------- 


Internal vertical motion for (dT/dx). < O and 


oe Pai WMS Tat Chesce sc lOno en = Seem Sk 


Numerical solution for (dT/dx),< O and ET 


cm Get Soe ee ee eee 


39 


40 


41 


42 


44 


45 


46 


48 


49 





ACKNOWLEDGEMENTS 


The author wishes to express her thanks to Dr. Robert L. Haney 
for his guidance given throughout the preparation of this thesis. 

Special appreciation is given to the night-shift personnel of 
the W. R. Church Computer Center and to a number of understanding 
friends who have been patient with me during the writing of this 


thesis. 





I. INTRODUCTION 


In recent years scientists have found that highly energetic, rotat- 
ing water masses exist in the ocean below the mixed layer. Eddies, as 
these masses have been called, appear to be similar to the high and low 
pressure systems in the atmosphere. Their origin and the extent of their 
influence on the temperature structure of the ocean is not fully known. 
Meteorologists as well as oceanographers are showing an increasing in- 
terest in these eddies as a means of obtaining a better understanding of 
the parameters of air-sea interactions. 

In 1970 a group of scientists began preparations that would lead to 
documentation of the existence of eddies. The Mid-Ocean Dynamics Experi- 
mene (MODE) was a multi-national experiment conducted in the area between 
Bermuda and Florida. Through extensive measurements of currents, tem- 
perature, salinity and pressure the scientists discovered the existence 
of two types of eddies; deep water eddies and highly energetic eddies in 
the upper 1500 meters of the ocean. 

It is the purpose of this paper to study the possible influence of 
the upper-layer eddies on the thermal structure of the ocean following 
the methods used in previous studies of the oceanic thermocline. These 
thermocline theories are all based on the assumptions of geostrophic and 
hydrostatic balance, incompressibility, neglect of salinity and some 
form of the steady state first law of thermodynamics which in most 


general form is 


oT oT oT O°T 
5 os Oy 6b ea aay (1) 





Table I gives a list of symbols and constants used in this study. The 
general balance expressed by (1) is a balance between horizontal and 
vertical advection and vertical diffusion of heat. Stommel and Webster 
(1962) have considered the simple barotropic case which balances verti- 
cal advection and vertical diffusion. Robinson and Stommel (1958) 
ignore the effects of east-west temperature advection (term one of (1)) 
while Welander (1971) considers a balance of all the terms in equation 
(1). In each of the above studies, the effect of mesoscale heating was 
ignored. 

The plan of this study is to consider the simplest model in which 
the effect of heating by mesoscale eddies given by the term, Swr), 
can be added to (1). If the energy source for the eddies is available 
potential energy, then w'T' and its vertical derivative will not be 
zero. Since the eddies have been verified by the MODE data, it there- 
fore becomes important to examine the effect that this term may have on 
solutions to the thermocline problem. For this study, the correlation 
term, owt") was specified according to theoretical results of Gill, 
Green and Simmons (1974) for a baroclinically unstable eddy of wave- 
length 322 km. 

Since the solution to the heat equation used for this study depends 
on the amplitude of w'T', its profile with depth and on the wind forced 
Ekman vertical velocity, cases were run to test each of these factors. 
Analytical and numerical solutions are found for the cases in which the 
Ekman wind and then the mesoscale heating are set equal to zero respec- 
tively. A numerical solution is found when both the wind forcing and 


the internal heating are included. In all cases horizontal temperature 


106 





TABLE I. Symbols and Values of Constants 


time 

vertical velocity 
perturbation temperature 
radius of Earth 
latitude, 36°N 

longitude 


angular speed of Earth's rotation, 27 day? 


- 1 


thermal expansion coefficient, 2.0 x hom Ron 
coriolis parameter, 282 sin ¢$ 

northward velocity component 

eastward velocity component 


: : ee =i. 
vertical eddy viscosity and conductivity, .5 om? sec 


iy 





advection is ignored. The results obtained in this study indicate 
that the effect of baroclinically unstable eddies on the ocean tempera- 


ture profile is significant. 


eZ 





II. EQUATION AND MODEL DEVELOPMENT 


A. HEAT EQUATION 


The basic heat budget equation for the ocean is 


ay 
Fie oy Oe (2) 


where Q represents the convergence of net incoming and outgoing radia- 
tion. It is assumed that radiative exchange takes place in the first 

few meters of the ocean. The model designed for this study is such that 
the level of z = 0 is at the bottom of the surface Ekman layer. This im- 


plies that within the ocean itself 


oll 
ae 0 (3) 


Expanding (3) and writing it in flux form gives 


= +V- (WW) =0 (4) 


Let V = V+ WW! and T=T+ TT!" where at this point the primes denote 

all turbulent motions including the baroclinic mesoscale eddies. A time 
average sufficient to remove all the turbulent motions (say 120 days) is 
performed with the result 


dm - = arr 
ae? Vr+V Cer) =o (5) 


3 





Expanding this into component parts and making the simplifying assump- 
tion that advection and heat transfer by turbulence in the north-south 


and east-west directions can be ignored gives the following equation 





is 


Ce 
me eee (6) 


ie 
trie 
+_ 
= I 
a? 


Zz 


The turbulent vertical transport of heat, (w'T'), can be thought of 
as being accomplished by two distinct physical processes. The first is 
through transport by baroclinic mesoscale eddies represented by the 
term (w'T') where the subscript of mesoscale will be assumed 

mesoscale 
for the rest of this paper. The second means of heat transfer is by all 
other turbulent motions. This transfer is assumed to be adequately re- 
presented by vertical diffusion of heat down the mean temperature 


gradient; «= - The final equation to be solved is therefore 


oT - OT eg pear oT -_ 
a Sage ge =. oe = 0 (7) 


where w'T' now represents the heat transport by mesoscale eddies. 


B. M EQUATION 

This study follows the approach of Welander (1971) who considered 
the solution to a similar heat equation for the thermocline region, 
but without the heating due to mesoscale eddies. The steady state, 
geostrophic, hydrostatic equations in a region away from side bound- 


aries are 


1 


Rcosd St Sees 


-fv =- 


Wir 


P, , PL = gar 


ll 
= cost {ud + (veos$) 4} +w = 0 


u 4 
———_ T, + —T, + = 
Rcosdé “A RO Wes Stee 
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where the variable P and T are perturbation pressure and perturbation 
temperature. Salinity is neglected and the assumptions of a spherical 
earth, ocean depth H small compared to the radius of the earth and con- 
stant density in the horizontal direction are made. Following Needler 
(1967), Welander derived a non-linear equation for pressure of fourth 
order and second degree. Noting the difficulty in solving the equation 
for pressure he defined a new function 
Z r 
M(A,$,2) = f pac + 2nReein$ woo (8) 
re) re) 
which is a first integral of pressure. The function also incorporates 
a boundary condition since Wa is the vertical velocity at the bottom of 
the surface Ekman layer. Welander points out that in general (u,v) # 0 
at the bottom thus implying a barotropic component. The M function pro- 


vides the following definitions: 


———_ y= eee" 
22Rsind oz 22Rsindcosd Az 
ail 7 2 a PP “e = 
22R sin > J 


Substituting these definitions into the heat equation (7) developed 


for consideration in this thesis, an M equation similar to Welander is 


obtained 
<a ) + a a M -KM + eck etm =e (9) 
ot 22Z ee A Z2z ZZZZ g Oz 
22R sin 
Letting 
Mo = 6 (10) 


As) 





and assuming 


K £ K(z) 
there results 
ie a —— 4,9, =U oes, + go wD") = 0 (11) 
2NR sin o> = 


The first term of (11) 1s proportional to the local time change of tem- 
perature. The second term represents the vertical advection of mean 
temperature while the third term is the vertical diffusion of the mean 
temperature. The last term is the heating due to mesoscale eddies. 
Numerical solutions to equations (10) and (11) for a steady state were 
obtained using the initial value approach. Analytic solutions were also 
obtained in the simplest cases. Before discussing analytic and numerical 
SOlutions to (10) and (11), a detailed description of the boundary condi- 


tions will be given. 


C. BOUNDARY CONDITIONS 

The physical model chosen to test the influence of ocean eddies is 
shown in Figure 1. The model is tested in a region of the Pacific Ocean 
away from any strong boundary currents, along 36°N. The domain was 
chosen to be large compared to the size of an eddy and to be sufficiently 
small so that the atmospheric forcing could be considered approximately 
uniform. A flat bottom was assumed for the two dimensional model. For 
initial test cases the boundary conditions were independent of longitude 
making this a one-dimensional model (solution independent of x). 

Since equation (9) is of fourth order, four boundary conditions can 


be applied to the problem. An advantage of using Welander's M function 
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‘Figure 1. Two dimensional model along 36°N with Az = 40m and 
Ax = 80m. 
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is that part of the boundary conditions are incorporated in the function 


itself. In practice it was convenient to solve (10) and (11) as the 





system 
Oe - oF ae) 
38 1 *) 
— + ———— M8 - KO + ga — (w'T') = 0 (a) 
dt 20R*sin*¢ Az ZZ OZ 


These are solved conveniently with boundary conditions on My at the 
ocean top and bottom in equation (12) and on 9 at the ocean top and 
bottom in equation (13). 

The model has a rigid lid with w prescribed at the top and bottom. 
ae a= 0, M, at the surface becomes 29R- sino Wo where a is a pre- 
Scribed vertical motion produced by the wind at the bottom of the 
surface Ekman layer. At z= -H , the boundary condition is w=0, 
which requires My = 0. The thermal conditions at the top and at the 
bottom are to prescribe the temperatures To and T respectively. 


=H 


Thus, at the surface 86 = ga T, and at the bottom, 9 = ga soe : 


D. EDDY HEAT TRANSFER 
In this thesis no attempt has been made to parameterize the eddy 

transfer of heat, w'T', in terms of the dependent variables. No direct 
measurements of the heat transfer were available so it was necessary to 
resort to theory to obtain an estimate. In a paper by Gill, Green and 
Simmons (1974) a theoretically determined rate of conversion of avail- 
able potential energy to eddy kinetic energy is developed as a function 
of z in a vertically continuous model. The values of w'T' presented by 


Gill et al are those that would exist if the mesoscale eddies were the 
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result of in situ baroclinic instability. As such, the w'T' which is 
used in this paper is based on a given horizontal and vertical profile 
of temperature. Figure 2 from Gill et al shows w'T' as a function of 
depth for eddies of four different wavelengths. The amplitudes asso- 
ciated with these profiles are based on the assumption that all of the 
available potential energy in the mean flow is converted into the 
eddies. This study can be thought of as a test of this assumption as 
Meeeedc a test of thé in situ baroclinic instabilaty theory. Figure 3 
Shows the mean velocity which determines the temperature profile. 

For analytic treatment in this thesis the heat flux due to meso- 


scale eddies was expressed as an exponential function. 





Woloe= Re (14) 


where A, .0357 °C cm Fare is the amplitude and D , 100m, is the 
characteristic depth of energy transfer. The profile of wit") 1s 

shown in Figure 4. This corresponds in amplitude and profile to that 
shown in (b) of Figure 2. The large amplitude indicates an extremely 


strong upward flux of heat in the upper 500 m of the ocean. 


Lg 
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Fig. 12. The rate of conversion of availadle potential eneryry to eddy energy, shown 2s a function 

of depth. The scale is chosen so that the integral with respect to depth (ig m) gives 311 Wo m7*9-. 

meeecasts shown arejil« -= a, £ = 0-033 kn [i]s = 01k = 0-005 hai? {2} < — 0-0235 km. 
(3}4 == 0-02758 kari. 


Figure 2. ‘Profile of the rate of conversion of available potential 


energy to kinetic energy. From Gill, Green and Simmons 
(1974). 
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Fig.7. Nfean velocity profiles U(z) for which a stability analysis has been carried out. Forclarity, 
profiles {4j and [$] have been displaced to the nghtand profile (I}is drawn twice. Note that the scale 


we 


Of sis timiced to fess nan half the sotal depth. The mean currents below 2000 m are very small. 


Figure 3. Mean velocity profiles from Gill, Green and Simmons 
(1974) used to determine temperature profiles. 
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IIIT. ANALYTIC AND NUMERICAL RESULTS 


The solutions to equations (12) and (13) are dependent on the inter- 
nal heating and the surface boundary conditions of temperature and verti- 
cal velocity. In order to aid in the interpretation of the more general 
numerical solutions, analytical solutions were obtained for a variety of 
the simpler cases first. 

The first case to be examined was that of no Ekman forcing and the 
next case was that of no internal mesoscale eddy heating. The final case 
was the most general case which included both Ekman forcing and heating 
by mesoscale eddies. Both upward and downward vertical velocities were 


tested where applicable. 


A. VERTICAL MOTION AT THE SURFACE EQUAL TO ZERO 
i. Analytical Solution 
In the first case the vertical motion at the bottom of the surface 
Ekman layer, Wor was set equal to zero. If the surface temperature, To! 
is independent of x, the solution, 6, is independent of x. Therefore by 
equation (12), My = 0 everywhere. The steady state form of equation (13) 


thus becomes (using (14)) 


2 er 
22-2) = 2 Q7/? (15) 
Oz 


This is the equation for a local, one-dimensional model without vertical 


advection. The analytic solution to (15) was obtained by integrating the 
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diffusion term twice and applying the boundary conditions of T(0) = oe 


27D 


meee y= ae) (16) 


Da) es Ge trad (le 

The solution (16) 1s shown in Figure 5. The temperature is 
forced to a value of -554°C at approximately 360 m below the surface. 
The extremely unrealistic nature of this solution indicates that if the 
open ocean contains mesoscale eddies with the above heating profile, 
then the thermal balance cannot be maintained by simple downward heat 
diffusion with kK = .5 com? sec '. 

2. Numerical Solution 


In order to test the numerical model which will be used later 


for the more general case, a numerical solution to 


— ar ce t ¢ 
ees eae ra ay oo 





was obtained. The numerical method is the initial value integration 


starting from an initial temperature field given by 


Tey) eee eee ae 


1 5 2: (18) 


The initial profile was defined to meet the stable condition of 


2s O. For this to hold the following constants were chosen; T) = 
meee C, T., =o Cc, ds = 900 m and d. = 100 m. Figure 6 shows this 
profile. 


The general numerical finite difference scheme to solve (17) 


used a centered differencing in space and time except for the diffusion 
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Analytic solution to the case 
equal to zero. 


25 


1000 


S00 


-200 -100 


in which wind forcing is 


nC 








_ _ OO Oo oe ae _ = — oa 


1000 


500 


-500 


-1000 


-1500 


-2000m 


Figure 6, Initial temperature (°C) profile for numerical model. 
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term which utilized a forward time step. The finite difference analogue 
is 


Te = ie + 2ht (kT * - go (wt) ) (19) 
where At = 10 days, Az = 40 m, Ax = 80 km. The final solution was 
reached after integrating for approximately 2000 years and is in agree- 
ment with the analytic solution. Intermediate solutions showing the 
approach to steady state are shown in Figure 7. It can be noted that 
the forcing term, (wT), immediately works to cool the temperature 
in the top 100 meters and then the diffusion term acts to redistribute 
the cooler water throughout the entire layer. Solutions to the baro- 
clinic case forced by a thermal gradient are considered in section D 


below. 


B. INTERNAL FORCING EQUAL TO ZERO 
1. Analytical Solution 
In the second case to be tested the heating due to mesoscale 
eddies, w'tT!, was set equal to zero. This required that there be a 
balance between heat diffusion and vertical temperature advection for 


steady state. From (12) and (13) the equations to be solved were 


20R*sin“¢ 
W = T 


ga. Zener so) 
2 

Kitevs (21) 
92° & 


If the surface conditions, To and Woe are independent of x, equation (2) 


reduces to Wee = 0 with solution 


z co 
w= wo(l + =) (22) 
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Figure 7. Intermediate numerical solution to case in which 
wind forcing equal to zero. a) 28 years b) 55 
years c) 82 years d) 165 years. 
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The equation to solve is thus 


2 
K2Fewa+y Ss (23) 


927 





with boundary conditions T(0) = To and T(-H) = Dy , Sirs 1S the local, 


one-dimensional equation with Ekman forced vertical motion of 10 cm ges 


but no heating due to the mesoscale eddies. 


Let y -= , then (23) becomes 


y' — _2 + Z 
. Z (1 2 
Wo 2° 
=— + —) + 
Thus, Qny r: (Z a Ci 
Ww 
Oo 2 
oT mene + Z /2H) 
ena y a ae = Cc. e 
Integrating from -H to Z; 
Ww 
Oo 2 
ain + T7422) 
T(Z) - Tha Bio e dn 
-H 


where mn is the variable of integration. Evaluating c. by requiring 


Z 
T(O) = z gives the final solution 
“ 2 
pone ren Une) 
T(Z) = Tha + rs w s e dn (24) 
f —(n+n /2H) -H 
K 
e dn 


eH 


The solution is not in closed form but is given by an integral 


Which must be calculated numerically. The results of equation (24) 
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being obtained by numerical integration using the trapezoidal rule are 
shown in Figure 8. Notice that the stratification is confined to the 
upper 500 m. 

Stommel and Webster (1962) have analyzed a similar case in 
which vertical heat diffusion and heat advection are balanced. Their 
solution in the barotropic case shows upward vertical velocity a linear 
function of depth, as in (22), but the temperature is isothermal. The 
Similarity transform used to solve their coupled equations does not 
allow vertical temperature stratification in the barotropic case. 

Figure 9 shows the case from Stommel and Webster (1962) in which the w 

and T lines denoted by 0° (implying a zero horizontal temperature gradi- 

ent) are for the barotropic case. The stratified temperature profile 

obtained in this study is allowed because equations (12) and (13) are 

solved without any restrictions imposed by a similarity transformation. 
2. Numerical Solutions 

In order to check the numerical model, the same case of no heat- 
ing by mesoscale eddies was solved numerically. By assuming that the 
local temperature change with time, =, is not equal to zero, an initial 
value time integration was run on the computer. The conditions are 
identical with those leading to equation (23) so the equation to inte- 


grate was 


2 
ee = re eee (25) 
dt 22 fe) H’ dz 
Z 
The finite difference analogue to this was 
— 4 —_ 
Tie ian eon ene aw (ler) (26) 
ZZ re) H Z 
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Figure 8. Analytic solution to the case in which heating by meso- 
scale eddies is equal to zero with La 10 cm day71l. 
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The model was run with a time step of ten days for 1000 years. A steady 
State was reached after 500 years for the case when we > O (upward 
motion). The numerical solution agrees with the analytically obtained 
solution shown in Figure 8. Since the numerically and analytically ob- 
tained solutions agree in these simple cases, it is reasonable to trust 
the accuracy of the numerical solution in the more general cases con- 
Sidered below. 

A numerical solution to the case with no internal heating was 
run with downward vertical velocity of the same magnitude as used before; 
LO cm oy Figure 10 shows an isothermal condition to great depth 
with a strong, deep thermocline. It appears that the vertical tempera- 
ture advection term works in the same sense as the vertical diffusion. 


Such a profile is clearly not observed in the real ocean. 


C. GENERAL CASE 

When both physical effects of heating by mesoscale eddies and wind 
forcing are included recall that the equations to solve are (12) and 
(13). In the case in which To and w, are independent of x, the solution 
1s independent of x, and so w is given by (22). Equation (13) becomes, 


with substitution of (22), 





30 Z ST 
Reo ms “oe * 8G, 


ll 
Oo 


(27) 


The steady state version of this equation is similar to (23) with the 
addition of the non-homogeneous term and an analytic solution is 
possible. Due to the complexity of the anticipated solution to equation 


(27) an analytic solution incorporating all three terms was not obtained. 
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Figure 10. Numerical solution to the case in which heating by meso- 
scale eddies is equal to zero with ie -10 cm day ~. 
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Only the initial value numerical method of solution was considered. 
Using a time step of ten days steady state solutions, shown in Figures 
11 and 12, were reached in approximately 500 years for the cases of 

Wo > O and Wo < 0 respectively. 

When Wo is upward (Figure 11) the initial temperature change due to 
the forcing function is similar to the case when ee QO (Figures 5 and 
7). ‘The first 250 m are cooled rapidly but then the influence of the 
vertical advection slows the cooling which is then confined to the upper 
region. The profile is more realistic than that given in Figure 5 for 
the case w = 0; however the temperature is still an unrealistic -174°C 
fe200 Mm. 

When Ww is downward, the result 1S Seen in Figure 12. Again cooling 
takes place quite rapidly due to the influence of the mesoscale eddies. 
The downward vertical velocity works with the vertical diffusion so the 
deeper waters also experience cooling. It has been shown that none of 
the local solutions which contain heating by mesoscale eddies are realis- 
tic. Further conclusions will be reached after examining the baroclinic 


cases below. 


D. BAROCLINIC CASES 
1. Internal Heating Equal to Zero 
The final series of test cases were designed to determine the in- 
fluence of heating by mesoscale eddies in a baroclinic ocean. This test- 
ing was accomplished by comparing the results obtained with and without 
the mesoscale eddy heat transfer term, (wt). The initial value 


approach to the solution was used with a temperature field having a five 


degree surface gradient. 
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Figure 11. Numerical solution to the general case with 
Wo = +10 cm day~1. 


216) 





1000 


S00 


-—250 = 200 -150 -100 -50 





Figure 12. Numerical solution to the general case with 
wo = -10 cm day71. 
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The equation, 


(Se (28) 


Z 
T(x,Z) = (T) # 5. (x/L))e 5 


where L is the east-west extent of the model (2000 km), was used to 
generate an initial field with either a positive or a negative east- 
west gradient. The values of T(x,0) and T(x,-H) obtained from (28) 
constituted the fixed thermal boundary conditions. With the addition 
of an east-west temperature gradient, the internal vertical motion is 
now a function of depth and temperature. Because equation (13) is now 
non-linear, My must be calculated at each time step. 
a. Positive Temperature Gradient 

Figure 13 shows the initial temperature field with a posi- 
tive gradient; warmer water to the east. The numerical solution to 
equations (12) and (13) is shown in Figure 14 for the case when the sur- 
face Ekman forcing is greater than zero, = 10 cm a An unrealis- 
tically deep isothermal layer is depicted. This is also the case for 
Ue less than zero. A strong interior downward vertical motion, shown in 
Figure 15, is created by the temperature gradient in both of the above 
cases. This sinking motion works with the diffusion to cause a packing 
of the isotherms near the bottom of the ocean. With no other forces act- 
ing to counter the double influence of the downward heat advection and 
diffusion, the resolution of the model is insufficient to handle the 
strong gradients that finally develop. 

b. Negative Temperature Gradient 
The initial field with a five degree negative east-west 


gradient at the surface is shown in Figure 16. When the model was run 
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Figure 13. Initial temperature field with a 5°C positive gradient 
in the east-west direction. 
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Figure 16. Initial temperature field with a 5°C negative gradient 
in the east-west direction. . 
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with cyclonic wind conditions generating an Ekman layer upward motion 
of ir 10 cm ge the resulting interior vertical motion is strongly 
positive throughout the model (Figure 17). The resulting vertical ad- 
vection is stronger than the downward diffusion and causes a packing of 
the isotherms near the surface. 

When anticyclonic wind conditions, generating downward ver- 
tical motion at the baSe of the surface Ekman layer, are applied the 
interior vertical motion becomes positive only below 250 m. Its maxi- 
Mum value is smaller than in the previous case in which w, was positive. 
Figure 18 shows a profile of the vertical motion taken from the central 
region of the model. For this case the vertical temperature advection 
does not dominate the downward diffusion and only a slight packing of 
isotherms occurs. A steady state solution, reached after 50 years, is 
Shown in Figure 19 with a well defined thermocline present. This is 
clearly representative of the typical situation below an atmospheric 
subtropical high pressure system with near surface temperatures gradual-~ 
ly decreasing eastward from the center of the oceanic anticyclonic gyre. 

Without internal heating the effect of a temperature 
gradient is as follows. When the east-west surface temperature gradient 
1S positive, the internal vertical motion becomes negative regardless of 
the surface value for Woe The resulting temperature profiles are similar 
to Figure 10. A negative surface gradient produces positive interior 
values for w and the baroclinic solutions without internal heating are 
Similar to Figure 8. 

2. Influence of Internal Heating 
The influence of heating by mesoscale eddies was tested for the 


baroclinic case with a negative horizontal temperature gradient and 
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Figure 18. Internal vertical motion for (0T/dx),< 0 and w_ = 
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Figure 19. Numerical solution for (dT/dx).< 0 andw< 0. 
Temperature is in °C. (d/dz) (w'T')= 0. 
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positive Ekman forcing. Equations (12) and (13) were solved using the 
initial value approach with the interior vertical motion calculated at 
each time step as a function of depth and temperature. After 40 years 
of integration, using a At of ten days, the result was a vertical motion 
profile shown in Figure 20. Recall that in the local solution case 
which included the eddy term with ie > 0 the interior value of w was 
simply a linear function of depth (Section III, C). 

The accompanying solution for temperature in the baroclinic case 
is shown in Figure 21. In this case the influence of the eddies was to 
act in the same sense as the interior vertical motion and cool the near 
surface layer to an unrealistic -142°C. This is comparable to the local 
case with ae > 0 (Figure 11), in which the minimum temperature reached 
was -174°C. It is evident that the baroclinicity acts to moderate the 
cooling influence of the mesoscale eddies but only to a small extent. 
Sweet al showed that the amplitude used for w'T' in this paper was re- 
Gumeea in order that in situ baroclinic instability produce eddies of 
the intensity observed in the open ocean. It is of interest to determine 
the sensitivity of the solution to the amplitude of w'T'. est cases 
were run in which the amplitude of w'T' was varied. It was found that 
reducing the amplitude by two orders of magnitude was required to produce 
a mean temperature profile close to that observed in the real ocean. 

Cases were also run with internal heating and the following con- 


oT oT oT 
ditions; 7 ~ < 0 withw <0, 7° > 0 with w > 0 and =~ > Owithw < 


ox Ox fe) Ox 
O. In all of these cases the resolution of the model was insufficient 
to handle the strong vertical temperature gradients which developed near 
the sea surface. Finer vertical resolution and thus smaller time steps 


would be required to obtain a numerical solution for any of the above 


three cases. 
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Figure 21. Numerical solution for (OT/Ox), < O and woe 
10 cm day71. 


49 





IV. CONCLUSIONS 


It is obvious that the real ocean does not freeze as a result of 
vertical heat transport by mesoscale eddies. There are a number of 
possible explanations for the excessive cooling observed in some of the 
test cases in this study. In deriving a model that was simple enough 
to test the influence of eddies, it was necessary to make numerous 
assumptions. With the motion independent of y and with the temperature 
Specified at the top and bottom of the model, the Reo seae eddy heat 
transport could only be distributed vertically by the term Zo wt). 
Were there an advective motion in the meridional direction, the north- 
ward eddy heat transport, vite, could have been distributed by the 
term air) such as to produce much less drastic effects and more 
realistic temperature profiles. It would be of interest to include 
horizontal temperature advection, v =, and aT) and then run the 
model testing the sensitivity of those terms. 

Since the mesoscale heating by transient eddies was determined 
from an analytic consideration of energy transfer rate by Gill et al 
(1974), it is possible that the amplitude of that transfer was overesti- 
mated. As stated before, Gill et al are estimating a maximum eddy 
strength by assuming that all of the mean flow available potential 
energy is converted into eddy energy. Since the earth's rotation acts 
aS a constraint to energy conversion and since the most favorable condi- 


tions for eddy growth are found where the isopycnals slope upward 


towards the equator, l1.e., colder water near the equator, a complete 
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conversion of mean available potential energy to eddy kinetic energy 
May not be possible. Using the amplitude for w'T' from Gill eteal, 

the unrealistic results from the model strongly suggest that the eddies 
Saeemvea in the open ocean are not the result of in situ baroclinic 
instability. 

The results of this thesis give a preliminary indication of the 
potential magnitude of influence that energetic mesoscale eddies have 
on the temperature structure of the ocean. It is evident that the 
heat equation in the ocean should include mesoscale eddy heating terms 


1£ accurate and complete solutions are to be obtained. 
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